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Abstract. We test the halo dominated mass model, re- 
cently proposed for high-luminosity galaxies by Courteau 
and Rix (1998), against rotation velocity data for a large 
sample oi L ^ L^, spirals. We find that this model does 
not reproduce the general kinematic properties required. 
The alternative disc dominated model, on the contrary, 
easily fits the data. 
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1. Introduction 

In a recent paper, Courteau & Rix (1998; hereafter CR98) 
have used the residuals of the well known relationships 
of spiral galaxies, luminosity vs. velocity and luminos- 
ity vs. radius, to estimate the disc-to-halo mass ratio 
P at R = 2.2 Rd^. Hereafter, V2.2 = V{2.2 Rd) and 
Ropt = 3.2Rd. The argument of CR98 relies on the idea 
that, for a given mass, more compact self-gravitating stel- 
lar discs have higher rotation speeds than less compact 

— 1/2 

ones (according to a Rjj scaling) . This effect should be 
observed in the rotation curves, if the disc component is 
dynamically important, and should disappear if the halo is 
dominant. More specifically, the effect can be measured by 
correlating the deviation 91ogV2.2 that an object has with 
respect to the mean velocity-luminosity relation, with the 
deviation dlogRo that the object has with respect to the 
mean radius-luminosity relation. For disc dominated sys- 
tems, one should expect 91ogV2.2/91ogi?£) = —0.5, while 
for halo dominated ones it is likely that d\ogV2.2/d\ogRD 
^ 0. In large samples of ~ spirals, CR98 find an ensem- 
ble value of d\ogV2.2/d\ogRD^ 0, which is interpreted as 
evidence that, on average, the disc contribution to V2.2 is 
far from maximal (see below) . The actual value of P is es- 
timated by computing the final halo mass distribution by 
adiabatically contracting (Blumenthal et al. 1986) a proto- 
disc inside an initial standard CDM halo (e.g., Navarro et 

Send offprint requests to: P. Salucci; e-mail:salucci(§sissa. it 
^ The radius where the rotation curve of a thin exponential 
disc of lenghtscale Rd peaks (Freeman 1970). 



al. 1996) into a final configuration with Rd — i kpc. So, 
CR98 find /3 = OMto.n for galaxies with 0.3< Lb/L^< 3 
(virtually independent of the initial DM distribution or 
the exact disc scale length or the presence of a [realistic] 
stellar bulge). 

Before discussing this result, let us comment on the 
procedure adopted by CR98. Three crucial assumptions 
have been made: 1 ) the residuals of both relationships are 
essentially due to differences in the mass structure be- 
tween each individual galaxy and the "average" galaxy; 
2) the residuals are not affected by discrepancies between 
redshift distancies and true distances, observational errors 
and selection bias; 3) the two relationships are intrinsically 
log- linear and/or do not have hidden parameters. The va- 
lidity of these assumptions is far from certain. Moreover, 
the method, devised to determine the average value of /3 
in a sample of galaxies, is poorly suited to detect sys- 
tematic and/or random variations of f3 within the sample 
itself. While these issues do not affect the main thrust of 
the CR98 idea, it is likely that they strongly modify the 
interpretation. 

Independently of the above issues (to be considered in 
a separate paper), the CR98 claim of /3 ~ 0.36 in high-L 
spirals is quite surprising. (Much less so if referred to low- 
L spirals, e.g. Persic & Salucci 1990.) Direct modeUing 
of rotation curves (RCs) of high-L galaxies, based on the 
maximum disc hypothesis (van Albada & Sancisi 1986; 
Sancisi & van Albada 1987), is able to successfully fit the 
RCs out to < 2 Rd with /3 ~ 1 (see: Kalnajs 1983; van Al- 
bada et al. 1985; Kent 1986; Begeman 1987; Broeils 1992; 
Broeils & Courteau 1997). Best-fitting multicomponent 
mass modelling, applied to extended high resolution RCs, 
has also found disc-dominated solutions for high-L galax- 
ies (e.g.: Broeils 1992; see references in section 4.) The RC 
profile methods, that are powerful diagnostics in quanti- 
fying the disc-to-total mass discrepancy, consistently find 
that ^ L, galaxies as disc dominated, i.e. /3 ~ 1 (Persic 
& Salucci 1988, 1990; Salucci & Frenk 1989; Casertano 



^ log(L«/L0) = 10.2; see Rhee (1996) for the B-magnitudes 
of the CR98 sample. We adopt Ho = 75 km s~^ Mpc"\ No 
result of this paper depends on this value. 
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Fig. 1. The halo velocity profiles of generic HD mass 
models with diflFerent values of /? < 0.5 which reproduce 
high-L RCs (dotted Unes). Also shown is the CR98 halo 
mass model (filled circles). The solid line is the PSS halo 
RC for high-L galaxies, plotted to explicit the difference 
between the HD and DD halo RCs. Velocities are in units 
of V2.2, radii are in units of Rd. 



& van Gorkom 1991; Salucci & Persic 1997). Finally, the 
swing amplifier constraints, coupled with stellar popula- 
tion arguments, seem to favour (near-)maximum discs (see 
Athanassoula et al. 1987 and Bosma 1998). 

Despite this, CR98 are not alone in their claim of 
halo-dominated (HD) high-L spirals. Mass models, built 
from the properties of the stellar velocity dispersion, find 
low values of /? (Bottema 1997); the swing amplifier con- 
straints may also imply subdominant discs (Fuchs et al. 
1998). Moreover, the local Galactic disc column density of 
(50 — 7O)M0 pc^^, estimated by studying the z-motions of 
local stars as a function of their height above the Galac- 
tic plane, is interpreted as evidence of a dominant halo 
(e.g., Kuijken & Gilmore 1991 and Kuijken 1995; but see 
Sackett 1997 and Dehnen & Binney 1997). 

On the cosmological side, the universal CDM den- 
sity profile (Navarro et al. 1996, 1997) features only a 
small central region where the halo rotation curve Vh{R) 
increases with radius, while it has Vh{R) ~ const for 
R > Rd that, combined with V{R) ~ const observed 
in ~ L* galaxies, implies a dominant halo. 

The case of HD mass models for high-L spirals is in 
(indirect) disagreement with our previous work (Persic & 
Salucci 1990; Persic, Salucci & Stel, 1996 [PSS]). How- 
ever, in view of the strong CR98 claim, it is worthwhile to 
investigate the issue further by means of specific tests. 

Notice that, for high-L galaxies, whose RCs are well 
described by V{x) oc — Ax~°' , with x = R/2.2 Rd, c ~ 
0, A ~ 0.1, a ~ 1, the CR98 mass model represents the 
whole class of HD models. In fact, a property of the above 
rotation curve is that any HD model (i.e., with (3 < 0.5) 
will have a halo velocity profile, Vh{R) = \J {V^ — V^^^^), 
very similar to the one adopted by CR98, independent of 
P and independent also of the precise values of c. A, a 



(see Fig. 1) in high-L objects. Then, a test performed on 
the CR98 model applies also to any HD mass model of 
high- luminosity spirals. 

In detail, in this paper we will devise two tests that, 
in the spirit of Salucci & Frenk (1989), single out two 
clear kinematical properties of RCs, and we will then check 
whether these are compatible with the CR98 model; failing 
that, we will iterate the tests on a disc dominated (DD) 
model. 

The plan is the following: in section 2 we determine 
the synthetic RC of high-L spirals; in sections 3 and 4 we 
perform two tests concerning the inner and outer rotation 
curves; in section 5 we discuss the results obtained and 
their implications. 

2. Rotation Curves of High-L Spirals 

In this section we coadd individual spiral RCs into a syn- 
thetic curve. This procedure, pioneered by Rubin et al. 
(1985) separately for Sa, Sb and Sc galaxies, was exten- 
sively used by PSS. In particular, for a sample of ~1000 
optical RCs of mostly late spirals, PSS have shown that: 
(1) within each given luminosity interval, the scatter of 
RCs around their mean is comparable with the scatter 
of data points around a given individual curve (i.e., the 
cosmic scatter is comparable with the error); and (2) the 
sequence of galaxy luminosities determines a sequence of 
RC shapes. [A Principal Component analysis of the same 
sample (Rhee 1996; see also Roscoe 1998) confirms the 
above results.] Moreover, at a given luminosity, there is 
practically no difference between the individual profiles of 
(high quality) RCs and the corresponding coadded profile 
(PSS). This evidence indicates that, as far as the gross fea- 
tures of their mass structure are concerned, spirals within 
a suitably small interval of luminosity can be considered 
as essentially alike: the coaddition procedure, therefore, 
eliminates most of the observational errors, extracting the 
relevant information stored in individual curves. 

CR98 use two samples of galaxies. One comprises over 
a hundred ~ L» Sb/Sc galaxies from Courteau (1992); 
the other comprises several hundred late spirals from the 
Mathewson et al. (1992; hereafter MFB) sample with 
Mi < -20.25. We then select from Persic & Salucci (1995; 
hereafter PS95) the RCs of MFB objects in the same lu- 
minosity range as CR98. From this subsample we will gen- 
erate a synthetic curve by coadding the individual curves. 
The procedure used is the following: we select galaxies 
with Mi < -20.25 and i < 85° (the RC of objects of 
higher inclination are excluded because, due to obscu- 
ration effects, the Ha position-velocity curves may not 
correspond to true rotation velocities: see Bosma et al. 
1992 and Bosma 1995). We then find 554 objects, with: 
< Mi >= -21.71, < Rd >= 3.6 kpc. | 



The average surface luminosity profile of the objects in the 
sample is well represented by a thin exponential disc (PSS). 
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Fig. 2. (a) Synthetic RC of ^ galaxies (filled circles), 
with superposed the halo dominated CR98 model (solid 
line), (h) The disc dominated mass model (solid line) su- 
perposed on the synthetic RC of high-i galaxies (filled 
circles). In both panels, velocities are in units of V2.2 and 
radii are in units of Ropu the bulge, disc and halo velocity 
components are denoted by dotted lines). 

Following PSS we coadd the (normalized) velocity 
data and obtain a raw synthetic RC with about 12000 
measurements, V{-^^; 0.3^ 3).^ We then form the 
smoothed synthetic RC of ^ galaxies, V{r), by averag- 
ing the velocity data over 24 radial bins of size 0.05 i?opt 
out to 1.2 Ropt (see Fig. 2). V{r) is the average RC 0/ (es- 
sentially) the same sample whose residuals have been stud- 
ied by CR98. 

3. Test I. Mass Modelling of Inner Rotation 
Curves 

In this section we check whether the innermost (i.e., for 
RS 2 R/o) RCs of ^ L^, spiral galaxies are halo dominated. 
The test is straigthforward because for R^ Ro the veloc- 
ity profile of a Freeman (1970) disc is completely different 
from that of a post-infall CDM halo. 

We first fit V{r) with the HD mass model which has al- 
lowed CR98 to model the residuals c?logV2.2 and dlogRo 
(see Fig. 7 of CR98) for the same sample that we have 
used to construct V{r). In detail, the model includes: (i) 
an exponential disc with Ru = 3.5 kpc; (ii) a stellar bulge 
with velocity profile peaking at Rm/Ropt — 0.045; (Hi) a 
post-infall CDM dark halo (see Navarro 1998). For this 

* In the following: r = R/ Ropt- 



mass distribution Vh{r) has an initial r^/^ profile and 
then a rapid convergence to a very shallow slope: set- 
ting Vh{r) oc r^^''\ then: fc(0.15) ~ 0.25 and k ~ 0.05 
for r > 0.7. 

The CR98 model is plotted in Fig. 2 a alongside its in- 
dividual components. It has one distinctive feature: about 
70% of the maximum velocity is reached at r = 0.10. By 
comparison, the exponential thin disc reaches the same 
fraction of the maximum velocity at a much larger radius, 
r = 0.21. Notice that the step-like shape of the halo domi- 
nated RC is not due the bulge, which reaches its maximum 
velocity at the much smaller radius r — 0.05, but just re- 
flects the shape of Vh{r) that, for r < 0.4, is steeper than 
Vdiscir). As a result, we realize that the model is incon- 
sistent with the data, by predicting, for < 0.4, a 
curvature in V{r) very different from the observed one. 

This failure does not occur only when testing the coad- 
ded RC, but also when we check individual RCs of high 
luminosity spirals that result severely inconsistent with 
the curve generated by the gravitational field of a domi- 
nant r-2_halo [see the cases of NGC 2179 and NGC 2775 
(Corsini et al. 1998), of NGC 5055 (Thornley & Mundy 
1997) and of NGC 5033 (Thean et al. 1997), and a few 
more cases in Begeman (1987) and Rhee (1996)]. 

Next, we fit V{r) with a DD model. This should not 
be considered as a further test, but rather as the neces- 
sary check that the coadded RC can be reproduced by a 
'natural' mass distribution, i.e. the one which follows the 
distribution of light. Applying DD models to samples with 
large luminosity ranges (as in the case of CR98) requires, 
however, the luminosity bins to be smaller. Then, in this 
case we analyse the RCs of MFB objects (see PS95) having 
Mi < -21.80. (The reader is referred to PSS for the anal- 
ogous results concerning lower-L galaxies.) The present 
subsample comprises 267 objects with global properties: 
< Mi >= -22.43, < logLg 10.52, < V2.2 >= 210 km 
s"^ and < Rd >= 4.3 kpc. 

The DD mass model we use includes: (i) sl thin expo- 
nential disc with Ru = 4.3 kpc; (ii) a bulge, with a Hern- 
quist velocity profile (practically indistinguishable from 
that used by CR98): 

^ ^ ^ Rm 1 + i?/ Rm ' ^ ^ 

where Rm = 0.63 kpc is the radius where the bulge veloc- 
ity peaks; (Hi) a dark halo with velocity profile: 

V?Ar) = V\l) (1 - /3) (1 + a^) , (2) 

where V{\) — 215 km s^^, a = 5 is the "velocity core ra- 
dius" (in units of 2.2Rd), and f3 = 0.85 is essentially the 
stellar mass fraction at 2.2 Ru (the bulge contributes only 
6% of the mass at this radius). Our adopted DD model, 
therefore, is not a best-fit model, but is chosen such as to 
only differ from the previously used HD model by (mainly) 
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the disc dominance: the bulge is kept virtually frozen (in 
normalized units), and the halo velocity profile is a sim- 
ple fit to the velocity residuals. The DD model easily re- 
produces the data for 0.1 within their rms errors (see 
Fig. 2 6): at very small radii the detailed modelling of the 
bulge is crucial but is beyond our scope here. 

Then we can conclude that, inside ~ 2i?£), the syn- 
thetic RC of high-L objects does not trace the dark halo 
mass distribution: in fact, it does precisely trace the mass 
distribution of the stellar disc. 

4. Test II. The Inner vs. Outer RC Slope 

A second test to check the ability of mass models to fit 
the RCs relies on the analysis of the outer RC profiles. 
The shapes of RCs can be essentially described by two 
quantities: the inner slope V = dlogV/d\ogR\ii^ ^ and the 
outer slope S = [Vi2Ropt) - y(i?opt)]/l^(^opt)Tsee PSS). 
Used in combination, S and V can help discriminating 
among apparently similar mass models. 

Let us consider, for example, two radically different 
models that lead to curves approximately flat at i? ~ Ropt ■ 
the "light disc -I- halo" model, and the "dominant 
disc -I- constant-density halo" model. In the former case, 
(iip)2 ^^d^ and (^)2 are both small in ab- 

solute value (< 0.1) but are opposite in sign, thus 

V ^ YkcJ_'^ + Yl±^ ^ 

1/2 Vd^sc dR y2 14 

implies also 5 ~ 0; moreover, since the disc contribution 
for R > Ropt is negligible, then V ~ i5 (specifically, the 
CR98 model has: V = -0.04, 5 = -0.04). In the second 
case, (^)2 and (^)2 are both a factor 

3-4 larger than in the previous case: they approximately 
compensate each other in the range \Ropt ~ Ropt, but 
further out the disc term takes over because in- 

Vdisa dR 

creases with radius while tt-^td decreases: then the dom- 

Vh dR 

inant disc model has —0.1^ 0.1 and 5 < V, with 
-0.2< 5< 0. 

We now check these predictions versus the observa- 
tions. From the literature we select, following PSS, 19 in- 
dividual RCs of high-L galaxies with a reliable profile out 
to 2 Ropt- These 19 curves are representative of the sample 
of 554 used for test I. (See Figs. 2 and 3 of PSS: the shape 
parameters V, 5 for the sample of individual RCs out of 
which these 19 are extracted, are consistent with the shape 
parameters of the coadded RCs built from essentially the 
same sample used in test I.) Using these curves, in Fig. 3 
we see that very few objects have a pair ^, V consistent 
with HD models (shaded area) and none lies in the region 
predicted by the CR98 mass model; moreover, ^ 90% of 

^ The "light disc -I- constant-density halo" and "dominant 
disc ~\- halo" models lead to curves that are respectively 
increasing and decreasing at Ropt- 
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Fig. 3. The distribution of high-L galaxies in the 
inner/outer-slope (V, S) plane. Data (filled circles) are 
compared with the CR98 model (the empty square). The 
shaded area represents the HD region. 

objects lay well inside the disc dominated region of the 
J-V plane. 

Finally, the constraint 6 < V < 0, predicted by the disc 
dominated model, is fully satisfied: the observed values 
are < S >= -0.08 ± 0.007 and < V >= -0.025 ± 0.008. 
Again, these values are inconsistent with the prediction of 
HD models: -0.04 < V < (5 + e, with e = 0(10"^). 

5. Discussion and conclusion 

Could high-L spirals be completely halo dominated? This 
can be effectively checked out by looking at their RCs: in 
fact, if the disc contribution to the RC is very small at 
all radii, the RC profile will essentially refiect that of the 
dark halo; otherwise, the stellar disc component will leave 
a recognizable imprint on the shape of the RC. 

By analysing the data for 554 inner RCs and 19 outer 
RCs of ~ L* spirals, we have uncovered two clear dy- 
namical imprints of a dominant disc. The first consists in 
a shallow extended rise in the RC between 0.3 Rd and 
1.3 Rd, that results coincident with that predicted by the 
Freeman (1970) disc; the second consists in a net decrease, 
(10 — 20) km s~^ (i.e. ~ 10%), of the circular velocity be- 
tween one and two optical radii, again as predicted by a 
dominant disc. These features, predicted well before sta- 
tistically relevant data were available (Salucci & Frenk 
1989), are inconsistent with the predictions of HD models 
(including the CR98 one). It is then obvious that mass 
models, featuring light discs and essentially constant ve- 
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locity halos, give a poor fit to the data, while models fea- 
turing massive discs and DM halos with solid-body-like 
profiles out to i? ~ Ropt match the data. 

On the basis of the results of two tests, one performed 
on the synthetic RC of high-L galaxies (seconded by lit- 
erature results on several individual RCs) and the other 
performed on about twenty individual RCs, we then reject 
the CR98 claim for halo dominated high-L galaxies. ^ 

We conclude by emphasizing the existence of a charac- 
teristic length-scale Rl{M{) such that, for radii R < Rl 
the baryonic matter (with total mass Mb) accounts for 
most (e.g., > 80%) of the gravitating mass inside R, and 
for R > Rl the dark component (whose contribution is 
minor at inner radii) rapidly becomes a major mass com- 
ponent. From PSS we can estimate: Rl ^ 2 Rd ijffrrM^)'' , 
with fc ~ 1/4. For high-L spirals, such "inner baryon dom- 
inance" extends to most of the disc region. 
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